CLOSED NODAL LINES AND INTERIOR HOT SPOTS OF THE 
SECOND EIGENFUNCTION OF THE LAPLACIAN ON 
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r^) ' Abstract. We build a one-parameter family of S— invariant metrics on the 

' unit disc with fixed total area for which the second eigenvalue of the Laplace 

^ ' ' operator in the case of both Neumann and Dirichlet boundary conditions is 

simple and has an eigenfunction with a closed nodal line. In the case of Neu- 
mann boundary conditions, we also prove that this eigenfunction attains its 
maximum at an interior point, and thus provide a counterexample to the hot 
spots conjecture on a simply connected surface. This is a consequence of the 
stronger result that within this family of metrics any given (finite) number 
r~| i of 5^— invariant eigenvalues can be made to be arbitrarily small, while the 

■^^ ' non— invariant spectrum becomes arbitrarily large. 



00 



O 
O 



1. Introduction 

JZ. • A conjecture of J. Rauch from 1974 states that the eigenfunction correspond- 

ly-s I ing to the second eigenvalue of the Laplace operator on a domain with Neumann 

boundary conditions attains its maximum and minimum on the boundary. As has 
0^ I been pointed out in Q, for instance, this is related to the location of the points of 

^^ ' maxima of solutions of the heat equation (hot spots) and it is basically equivalent to 

saying that, for most initial conditions, these hot spots move towards the boundary 
as time goes to infinity. 

Recently, Burdzy and Werner [g gave a counterexample to this conjecture on a 

C^ ' domain in R^ with two holes and posed the question of whether it would still be 

possible to find a counterexample on a doubly connected domain, or whether the 

conjecture would hold in that case ~ note that it is not known if the conjecture 

holds on simply connected domains. 

^^ ' On the positive side, Kawohl pSJ has shown that for the case of domains of the 

H I form D X (a, 6), where D C R'^'^has a C°'^ boundary, the values attained by a 

second eigenfunction on the domain are less than or equal to the values it takes 

on the boundary. More recently, Baiiuelos and Burzdy proved that the conjecture 

holds in the case of some special domains which include, among others, convex 

domains with a line of symmetry Bj . 

This problem is related (although not necessarily equivalent) to the nonexis- 
tence of closed nodal lines of the second eigenfunction, where the nodal set of an 
eigenfunction is defined to be the closure of the subset of the domain where the 
eigenfunction vanishes. With the exception of Courant's nodal domain theorem 
which states that the nodal set of a fc"^ eigenfunction divides the domain into at 
most k subregions H , very little is known about the general structure of such sets - 
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see [0, H ^ . Courant's result implies that a second eigenfunction always divides the 
domain into exactly two subregions and a conjecture of Payne from 1967 for the case 
of Dirichlet boundary conditions states that any such eigenfunction cannot have a 
closed nodal line pq |. This was proved by Payne in 1973 under some symmetry and 
convexity assumptions on the domain [^. Within the last ten years there have 
been several developments regarding the existence or not of closed nodal lines of the 
second eigenfunction for the Dirichlet Laplacian - see M R nfl, O, U2, 14 , isl 




Although it has been shown that the conjecture holds for convex domains 
again it is not known whether or not there exist simply connected sets where the 
second eigenfunction will have a closed nodal line. So far the only known counterex- 
ample requires the boundary of the domain to have at least three components |10| . 
Note that, in the case of the Neumann problem, it is not too difficult to prove that 
on simply connected domains the second eigenfunction cannot have a closed nodal 
line|,|l|,|l|. 

The question of whether or not results similar to these hold in the case of man- 
ifolds with boundary has been raised by S.T. Yau - see, for instance, Problem 45 
in the Chapter Open problems in differential geometry in |21[| . Note that in this 
case it is not difficult to find a counterexample to the closed nodal line conjecture 
if one considers doubly connected surfaces. For this, it is sufficient to think of long 
cylinders, for which the nodal line is a circle located halfway between the top and 
the bottom. Incidentally, this also justifies the statement that the hot spots and 
the closed nodal line problems are not necessarily equivalent, as in this case the 
maximum and minimum are still attained at the boundary. As far as we know, the 
only results for the case of surfaces are given in the book by Bandle ^ , where it 
is shown that in the Neumann case and under some conditions on the metric the 
second eigenvalue is not simple and the maximum and minimum are attained at 
the boundary. It is stated in that book that one cannot expect for this to be true 
in general, but no counterexample is given. 

The main purpose of this note is to present examples of simply-connected sur- 
faces with boundary for which both the closed nodal line and the hot spots conjec- 
tures fail. In other words, for these surfaces the eigenfunctions associated with the 
second eigenvalue of the Laplace operator in both the Neumann and the Dirichlet 
case have closed nodal lines which do not touch the boundary. Furthermore, in the 
Neumann case, this eigenfunction can be chosen in such a way that it attains its 
maximum at an interior point. This shows that simple connectivity by itself is not 
a sufficient condition for the above conjectures to hold. 

These counterexamples are corollaries to the stronger result which asserts the 
existence of a family of S'-'^-invariant metrics with constant area on the unit disc for 
which the first m eigenvalues are simple, for any given positive integer m. For S"^- 
invariant metrics, the spectrum can be divided into an invariant part corresponding 
to the spectrum of the Laplacian acting on S'^-invariant functions, and a non- 
invariant part. The counterexamples given here are then made possible by the 
fact that we are able to choose a family of metrics for which the invariant and 
the non-invariant parts of the spectrum may be separated. More precisely, for any 
positive integer m we have that the first m invariant eigenvalues can be made to be 
arbitrarily small, while the first non-invariant eigenvalue becomes arbitrarily large. 
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2. Main result 

The main result from which the counterexamples are then a straightforward 
consequence is the following 

Theorem 1. There exists a one-parameter family A4 of smooth S^ — invariant met- 
rics on the unit disc with positive curvature and total area n, such that given any 
positive integer m and real number e there exists a subset M.^ of M with e on an 
open interval for which the first m S^ — invariant eigenvalues of the Laplace opera- 
tor both with Neumann and Dirichlet boundary conditions are smaller than e. On 
the other hand, the non-invariant spectrum remains uniformly bounded away from 
zero in Ai and becomes arbitrarily large in M^ as e goes to zero. 

Since the Gaussian curvature is positive, these surfaces can actually be isomet- 
rically embedded in R^ [1^. As we shall see in Section ||, the curvature of such 
metrics can be made to be arbitrarily close to zero except at the centre of the disc. 

From the proof of the theorem, it follows that the first m invariant eigenval- 
ues must be simple, and, in the case of the second Neumann eigenvalue, that the 
corresponding (invariant) eigenfunction is strictly monotone along radial lines and 
changes sign. We thus obtain the following 

Corollary 2. Given any positive integer m there exists a family M of S^ — invariant 
metrics on the unit disc with positive curvature and total area tt, for which the first 
m eigenvalues of the Laplace operator with both Neumann and Dirichlet boundary 
conditions are simple. In both cases, the eigenf unctions corresponding to the j^^ 
eigenvalue (j = 2,... ,m) are also invariant by S^ and have j — 1 nodal lines 
which are closed disjoint circumferences dividing the disc into j nodal domains. In 
the Neumann case the second eigenfunction may be chosen in .such a way that its 
(.strict) maximum is attained at the origin. 

The estimates obtained in the proof allow us to write a more quantitative version 
of this result which we mention in Section H. 

Counterexamples to this type of conjectures have usually been obtained by 
exhibiting a domain for which the second eigenvalue can be proven to be sim- 
ple [^, |l^, |l^ . If the domain has some symmetry, it will then be inherited by the 
eigenfunction and hence also by the nodal set. In order to prove the above results, 
we shall use a variation of this technique and consider one-parameter families of 
metrics on the disc for which we are able to prove that the m"^ invariant eigenvalue 
is smaller than the first non-invariant eigenvalue. 

A fundamental (standard) ingredient in the proof is the reduction of the original 
two-dimensional problem to a sequence of one-dimensional problems, which we 
achieve by means of standard polar coordinates. Since we are interested in estimates 
for higher eigenvalues, we then perform a change of variables in order to avoid 
function weights in the orthogonality conditions which appear in the corresponding 
variational formulations - see Section ^. At this point we should remark that 
it is also possible to proceed by means of a different technique using symplectic 
coordinates. An example of this can be found in Q, where the behaviour of the 
invariant spectrum for 5^— invariant metrics on S^ was studied. We shall briefly 
indicate in Section || how these two different coordinate systems are related in this 
case. 

Regarding the choice of metrics, we point out that with the proper parame- 
trization in isothermic coordinates (see Section g) , the eigenvalue problem on the 
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surface becomes equivalent to that on a flat disc of inhomogeneous density, and so 
our results also apply in that case. On the other hand, this suggests that the we 
build the family A4 by making the density much higher close to the centre than 
near the boundary, so that there is a strong resistance to the movement of the hot 
spots towards the outside regions of the disc. We remark that in the opposite case 
where the (radially symmetric) density increases as we move away from the centre, 
it is known that the hot spots conjecture holds 0. 



3. Abstract surfaces 

In this section we collect the main facts about abstract surfaces that will be 
needed in the paper. We shall follow the exposition in the book by Bandle Q very 
closely, and begin by introducing the concept of an abstract surface S. We then 
derive the expressions that will be used in the sequel, which include the Laplace- 
Bcltrami operator in S given in conformal coordinates and the variational formu- 
lation for the problem. We shall also indicate the expression for the Gaussian 
curvature of S. 

Let 13 be a domain in the {x, y)— parameter plane, and da^ the Riemannian 
metric in D defined by the quadratic form 



da-" 



E{x, y)dx'^ + 2F{x, y)dxdy + G{x, y)dy'^ 



Definition 3.1. A domain D C M^ with the Riemannian metric da is called an 
abstract surface and will be denoted by 5 = {D, da). S is said to be in its isothermic 
(or conformal) representation if da^ = p{x,y)ds^ {E = G ~ p and F =^ 0), where 
ds denotes the linear element of the Euclidean plane. 



(3.1) 



The eigenvalue problems that we are interested in are thus 

Agu + ju = in D, 
= or u = on dD, 



du 
'dv 



where A5 denotes the Laplace-Beltrami operator on S and v the conormal deriv- 
ative. We shall denote the eigenvalues of the Neumann and Dirichlet problems by 
pLj and Aj, j — 1,..., respectively, and always assume that they are written in 
increasing order. 

The Laplace-Beltrami operator is now given by 



1 


d 

dx 


dx dy 
W 
\ ) 


d 

dy 


dy dx 
W 
\ J . 



where W = \JEG — F"^ . In the case where S is in its isothermic representation, 
this expression simplifies to 



Ac = 



1 

p{x,y)' 
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where A now denotes the usual Laplacian operator in M~. This means that the 



eigenvalue problem (3.1) on 5 becomes 



/ -*- , Am + 7M = in D, 
(3.2) P^^ 

1^ = or w = on dD, 

As mentioned in the Introduction, this is equivalent to the problem of an inhomo- 
geneous membrane whose density is given by the function p. 

Finally, we note that in this case, the area of S and the Gaussian curvature are 
given by 

Ap{S) = / p{x,y)dxdy and X = -— A [log(p)] , 
Jd 2p 

respectively. The expression for the total curvature is then 
^{S) = -- / A [\og{p{x,y))] dxdy. 

From this point on we shall concentrate on the Dirichlet problem, and mention 
the necessary changes for the Neumann c ase in Section Q. 

We consider the eigenvalue problem (|3.2| ) in the case where the function p is 
radially symmetric, that is, p — p{r). In polar coordinates and after separation 
of variables, the eigenvalue problem (|3.2| ) then reduces to the sequence of one- 
dimensional eigenvalue problems 



(3.3) {r^y 



Xrpir) 

r 



(^-0, re (0,1), 



where fc^, k = 0,1,..., are the eigenvalues of the Laplacian on the circle. For 
fc = we have Neumann boundary conditions at and Dirichlet at 1, that is 
</3'(0) = </?(l) = 0. This corresponds to the invariant spectrum of the original 
problem, and the associated eigenfunctions are the radially symmetric functions Lp. 
Note that the i^"^ eigenfunction of the one-dimensional problem has i — 1 zeros on 
(0, 1), and thus the corresponding eigenfunction of the original problem divides the 
disc into i nodal domains. 

For positive values of k we have Dirichlet boundary conditions at both ends of 
the interval, and this now gives the non-invariant part of the spectrum which is 
made up of eigenvalues with multiplicity two. In this case, two linearly independent 
eigenfunctions are given by i^{r) cos{k9) and (p{r) sm(k9). 

To obtain the variational formulation corresponding to the eigenvalue problems 
above we consider the space C^(0, 1) with the inner product defined by 

{(p,ip) ~ / rp{r)ip{r)ijj{r) + rip' {r)ip' {r)dr. 
Jo 

Let now Hq be the Sobolev space which is obtained as the closure ofC^{0, 1) with 
respect to the norm induced by the above inner product. Then the eigenvalues of 



the spectral problem (3.3) when fc = are given by 



r[(p'{r)] dr 

(3.4) A° = A°(p)= inf -^f , j = l, 

rp(r)(p {r)dr 



v^"L I t \ 2, 
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where Hq j^ is the closure with respect to the above norm of the space of C°° functions 
on [0,1) with compact support on this interval, 

Hi, = ii^oj^i ^ {"^ • / rp{r)^,^^{rMr)dr - o| , J = 2, . . . , 

and ipj is the eigenfunction corresponding to A*^. 

In the case of the eigenvalues for each of the remaining problems we have 

(3.5) 

X^ = X^ip)= inf ^-^ -^ ^^2 : , j = l,..., fc = l,..., 



rp{r)ip {r)dr 

and where now the spaces iJg ■ are defined in a similar way as above but starting 
with Hq. 

As we are interested in higher eigenvalues, we have to take into account the 
orthogonality conditions which appear in the definition of the spaces Hq ■ . Since 
these conditions become much simpler to handle if the inner product considered 
has a constant weight function instead of rp{r,S), in order to show that X%^{p) is 
arbitrarily small within a family of metrics it is convenient to make a change of 
variables in the variational formulation of the problem. Thus, we want to find a 
function r = r(z) and a constant c such that 

r{z)p {r{z)) r'(z) — c 

and the interval (0, 1) is mapped onto (0, 1). This will be the case if we take, for 
instance, 

(3.6) z(r) = — / sp{s)ds and c — sp{s)ds — ^ 



cJo Jo 27r 



The Raleigh quotients in (3^) and (3_^) then become 

r-l 

-2 



r^z)p{r{z))[^P\z)]^dz / [r^{z)p {r{z))] ^ [i>{z)]^dz 






(3.7) -^—^ 1 + fc'^ 

^P^^' [ ^^{z)dz I ^\z)dz 



"'0 



for k = 0, . . . . This transformation fixes both endpoints of the interval, and the 
inner product is changed accordingly to 



1 47r2 '■I 



i4;^)=Ap{S) ^^dz + —j- r\z)p{r{z))^'4>'dz. 

We thus have that the orthogonality conditions used to define the spaces appearing 
in the variational problem are now formed with respect to the usual inner product. 
We shall use the same notation for these spaces as above, as it will be clear from 
the context which inner product is being used. 

We remark that it is also possible to choose the change of variables above in 
such a way that the problem becomes equivalent to that obtained in ||I[ using 
symplectic coordinates. We thus obtain a relation between symplectic coordinates 
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and isothermic coordinates for this particular case. To do this, we need the relation 
between r and z to be such that 

(3.8) r p{r) — ——— and r'(z) = rg{z). 

Here g : (—1, 1) -^ (0, +oo) is a function of the form 

1 

where h is a. C°°[—l, 1] function such that g > 0, and it has been assumed that the 
area of the full surface is An. We then obtain that the change of variables defined 

by 



g{z) ^ ^ + h{z) 

I yZ 



h{t)dt 



^i^)^\l^r—:^ ■''' \l-, — z ^>0' 



transforms (3.4) into 



l-e 



^ [ij'{z)]'^dz 



9{z) 



2/ 



Tp {z)dz 
The function p is then given from (B^) by 

p{r) 



although it should be pointed out that in general it will not be possible to obtain 
an explicit expression for it in terms of r. Note that because polar coordinates 
have a singularity at the origin, while symplectic coordinates have two singularities 
at both the North and South poles, we have to restrict o urse lves to the interval 
(— 1, 1 — e). In this case, the change of variables given by ( |3.6|) has the advantage 
that it keeps the interval fixed. 

As an example, consider the case of the standard sphere for which we have that 
h is the zero function and then 

4£(2-e) 



as expected. 



P(^) = 9 

[e+(2-s)r^]^ 



4. Proof of Theorem 1 



The first m eigenvalues of the original problem will be simple if and only if we 
have that Aj^ is smaller than A^' for all positive k. This follows from the fact that the 
invariant spectrum is the spectrum of a second order ordinary differential operator 



of the type in (3.3) with k = Q and thus all its eigenvalues are simple. Clearly it is 
enough to ensure the above condition for fc = 1, since \\ < Xi~^^ for all k. To do 
this, we consider the one-parameter family of smooth nrctrics corresponding to 

a 1 



p{r, S) 



r 



2 



log 
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for positive 6. The parameter a is included as a normalizing factor so that all 
elements in this family of metrics have fixed area equal to tt. 
Using the expression given in Section H we easily obtain 

and thus the curvature is positive and converges to on (0, 1] as 6 goes to 0. On 
the other hand, the total curvature is given by 

'^(^^ = TTs- 



Using the change of variables given by (3.6) we obtain 



r = '^6{e^/°' - 1), 
and the variational problem (3.7) becomes 



(4.1) 



A^- = Af (p) = inf 



+ ■ 



4a ^ 



l-e-'-lA [tP'{z)Ydz 



\ 



-0 {z)dz 



1 



a 



1 - e-^'' 



if) (z)dz 



[ tlj'^{z)dz 
Jo 



J-1, 



fc==0, . 



/ 



We begin by noting that since a is positive the second term is always greater 
than or equal to k^ /a so that we have 

A}(P) > -, 
a 

which becomes unbounded as S goes to zero. 

We shall now show that any finite number of invariant eigenvalues can be made 

to be arbitrarily small. We have that 

"1 , V /"i 



1 



'Z I OL 



\i)'{z)Ydz 



\i)'{z)fdz 



< 

ip-{z)dz f ^p^{z)dz 

)o Jo 

and so by the monotonicity principle derived from Poicare's principle it follows that 

A0(p)<a(2j--l)V, j = l,..., 
and thus, given a fixed eigenvalue A'^(p), 

lim A°(p) =0. 
The estimates obtained above allow us to conclude that if 

^ < g(2j-i)^_-^' i = I'--- ' 

then the j*'' invariant eigenvalue is smaller than the first non-invariant eigenvalue. 
This is a very rough estimate, since the logarithmic term is taking into account 
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the asymptotic behaviour of the eigenvalues, while the crossing between the two 
eigenvalues as the parameter i5 decreases must occur at a point larger than 1. Similar 
results can be obtained for higher eigenvalues. 

Due to the fact that the Dirichlet and Neumann eigenvalues satisfy fij < Xj for 
all integer j, the proof for the Neumann problem follows from the Dirichlet case. It 
is also possible to obtain some independent estimates by following essentially along 
the same lines as in the Dirichlet case, with only some minor changes. These have 
to do with the spaces considered, which should now be 

Hl^H^ and Hj = Hj_-^^n lip : J^ ipj-i{z)ip{z)dz ^ o\ 



in the case of invariant eigenvalues, and Hq in ( |3.5[ ) is now the closure of the space 
of C°° functions with compact support on (0, 1], with respect to the same norm as 
before. We may then proceed in a similar fashion to obtain the following estimate 
for the j^^ invariant Neumann eigenvalue: 

Ai°(p)<4a(j-l)V. 



From ( |3.3| ) with fc = we see that any solution ip of that equation with precisely 
one zero on the interval (0,1) and with ip{0) positive will be strictly decreasing on 
this interval. Hence its maximum is attained at the origin. 
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